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An ultra-cold Bose gas of two-level atoms can be thought of as a spin-1/2 Bose gas. It supports 
spin-wave collective modes due to the exchange mean field. Such collective spin oscillations have 
been observed in recent experiments at JILA with 87 Rb atoms confined in a harmonic trap. We 
present a theory of the spin-wave collective modes based on the moment method for trapped gases. 
In the collisionless and hydrodynamic limits, we derive analytic expressions for the frequencies and 
damping rates of modes with dipole and quadrupole symmetry. We find that the frequency for a 
given mode is given by a temperature independent function of the peak density n, and falls off as 
1/n. We also find that, to a very good approximation, excitations in the radial and axial directions 
are decoupled. We compare our model to the numerical integration of a one dimensional version 
of the kinetic equation and find very good qualitative agreement. The damping rates, however, 
show the largest deviation for intermediate densities, where one expects Landau damping - which 
is unaccounted for in our moment approach - to play a significant role. 



PACS numbers: 05.30.Jp, 32.80.Pj 



I. INTRODUCTION 



Collective spin oscillations are a general consequence 
of the quantum exchange between identical particles in a 
system where a macroscopic symmetry breaking exists in 
spin space From a condensed matter perspective, spin 
waves are most familiar in strongly interacting degener- 
ate Fermi systems, such as a ferromagnet. It is somewhat 
counter-intuitive, although well established both exper- 
imentally |§, |, f§ and theoretically §, §, @, |, §, that 
collective spin behavior can also occur in nondegenerate 
dilute spin polarized gases when the thermal de Broglie 
wavelength exceeds the effective range of interaction be- 
tween two colliding atoms. In these systems, a trans- 
verse spin wave is excited by applying an inhomogeneous 
magnetic field followed by a small tipping pulse. It is re- 
markable that the mean field generates collective spin dy- 
namics, but has no discernible effect on thermodynamic 
equilibrium properties since gn/ksT <C 1, where g is the 
binary interaction parameter, n is the density and T the 
temperature. 

Recent experiments at JILA jlO|, [ll] on a trapped 
87 Rb gas have revived interest in spin waves in dilute 
gases [|l2| |l^, [TI], These new experiments offer sev- 
eral interesting new features compared to the earlier ex- 
periments in spin-polarized hydrogen ^ 0], which take 
advantage of the technological advances made over the 
last twenty years in the measurement and control of cold 
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atomic gases. A prominent feature of the new gener- 
ation of experiments is the ability to take spatially re- 
solved measurements of the gas sample using absorp- 
tion imaging techniques. Another exciting advancement 
is the ability to cool the sample into the quantum de- 
generate regime, which permits the study of spin waves 
in a Bose-condensed gas at finite temperatures; this 
regime has never been investigated experimentally and 
has only received minor attention in the theoretical lit- 
erature 1 16, |l7, pH| . In this paper, however, we focus 
on the noncondensed regime relevant to the recent JILA 
experiments 0, |ll| , where the temperatures are approx- 
imately twice that needed for Bose-Einstein condensation 
Tbeo 

The JILA system consists of a dilute gas of 87 Rb atoms 
that have been optically pumped into the \F = 1, Mp = 
— 1) = |1) hyperfine state and are confined in a magnetic 
harmonic trapping potential. By applying microwave and 
radio-frequency radiation that couples to the |2, 1) = |2) 
state, atoms in the gas can be uniformly prepared in an 
arbitrary superposition of the |1) and |2) states. This 
system can be thought of as a spin-I/2 system by taking 
|1) as the spin-up state and |2) as the spin-down state. 
In Figure 1 we illustrate the corresponding Bloch vec- 
tor spin describing the internal state of the atoms. Note 
that because the magnetic field direction varies in the 
trap [ |l9| , the spin axis shown in Figure 1 is not isomor- 
phic with the coordinate axis describing the position of 
an atom, in contrast to the situation in spin-polarized 
hydrogen. 

In the JILA experiment, an initial 7r/2 pulse is ap- 
plied to tip the spins into the transverse direction v. The 
spin vector then precesses about the longitudinal w-axis 
at a rate proportional to the energy difference between 
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to a spin kinetic equation, and derive explicit analyti- 
cal expressions for frequency and damping of dipole and 
quadrupole modes in weak and strong coupling limits. 
We also numerically solve a one-dimensional model of the 
kinetic equation, and compare with the moment results. 



II. SPIN KINETIC EQUATIONS 

The Hamiltonian describing a single, trapped, two- 
level atom of mass m is: 



FIG. 1: Schematic diagram of the Bloch vector spin S. The 
spins of the atoms are initially pointing up, 5* = {0, 0, S w }, 
corresponding to all the atoms being in the state |1). After 
the 7r/2 pulse, the spin vector is pointing along the v axis, 
S = {0, S v , 0}, corresponding to the state (|1) + i\2})/V2. 



H = 



- — V 2 + C/ cx t(r) 
2 m 



l + -n(r)-f. (1) 



The first term in Eq. (Q) is the center of mass Hamil- 
tonian containing the kinetic energy and the external 
parabolic trap potential 
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hyperfine states. Due to the mean field and differen- 
tial Zeeman effects, the local frequency splitting between 
hyperfine states varies approximately quadratically with 
position. This inhomogeneity initiates collective spin dy- 
namics through the exchange mean field, the initial on- 
set of which gives the striking appearance of spin seg- 



regation |L0|, 12, |13|, |14J. The inhomogeneous frequency 
splitting can be made arbitrarily small to study the linear 
response of the system. Recently this technique was used 
to probe intrinsic collective spin oscillations |fj"l| . Stud- 
ies of such spin-wave collective modes provide us clearer 
physical understanding of the important role of the ex- 
change interactions in a dilute Bose gas. 

The experiments described in Refs. jnj |ll[ present 
spatially resolved images of spin dynamics in a gas. The 
density profile of either state is measured using absorp- 
tion imaging. Together with the Ramsey fringe tech- 
nique, integrated spatial profiles of the longitudinal spin 
S w and transverse phase <j) can be extracted from experi- 
mental data, as shown in the stunning images in Ref. [ pd| . 
This is in sharp contrast to the earlier hydrogen experi- 
ments, where a pulsed NMR technique is used to obtain 
the frequency of spin oscillation integrated over the entire 
sample. In this paper we show a theoretical prediction 
of the spatial structure of the spin dynamics that quali- 
tatively agrees with experimental observation. 

In this paper, we present the theory of spin waves in 
a trapped Bose gas using the moment method [^o[ fH| . 
which was originally developed to study collective den- 
sity oscillations in a trapped classical gas. The moment 
approach has also been applied to a rotating gas |||] 
and was recently generalized to treat a Bose-condensed 
gas at finite temperature to study the scissors modes 
p3[ . An advantage of this technique is that the solu- 
tion maps smoothly between the collisionless and hydro- 
dynamic regimes p0[. We apply the moment method 



where a — uj r /uj z . This part of the Hamiltonian is un- 
coupled from the internal, pseudo-spin, degree of free- 
dom, which is governed by the second term: f2(r) ' t — 
Qu(r)T u + Qv(y)t v +ri M )(r)f tu , where is a Pauli matrix. 
In the absence of an external coupling field, £l u = fl v = 
and fl w = A(r) [ll| is the frequency splitting between 
the two states (we go to a rotating frame to eliminate 
the large hyperfine splitting Whf frequency). We model 
binary interactions between particles by a delta-function 
pscudo potential describing elastic, spin preserving col- 
lisions, the strength of which depends on the hyperfine 



states V~ij(r,r') = gijS(r — r'), where 



with fly being the scattering length for collisions between 
atoms of species i and j. For 87 Rb, an = 100. 9ao, 
ai2 = 98.2ao, 022 = 95.6ao, where ao is the Bohr ra- 
dius [fiof . In the rest of this paper, however, we make 
the simplification that an = 022 = &12 = a, which is a 
reasonable approximation for 87 Rb. 

Several groups have previously worked out the funda- 
mental kinetic theory of a noncondensed dilute Bose gas 
with internal degrees of freedom, to describe spin waves 
in spin-polarized atomic hydrogen above Tbec |j- @H- [1 • 
Using a semiclassical approximation to describe atomic 
motion in terms of a phase-space distribution function, 
we obtain coupled Boltzmann equations for the atomic 
/(r, p, t) and spin cr(r, p, t) distribution functions: 



dt m pJ 2 p dt 



coll 

(3) 



at m 1 at coll 

(4) 

Equation (ra) has an implicit sum over the repeated in- 
dex i = it, v, w in the fourth term. The total density and 
spin density are obtained from the distribution functions 
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as n(r,t) = m(r,t) + n 2 (r,t) = J dpf(r,p,t)/(2irh) 3 

and S(r,t) — J dp<?(r, p, t)/(2irh) 3 respectively. Here 
the longitudinal component of the spin represents the 
relative density S w = ti\ — ri 2 and the transverse compo- 
nents S u and S v describe the real and imaginary parts 
of the internal coherence. The center of mass effective 
potential is U n (r, t) = f/ e xt( r ) + 3.gri(r, t) /2. The effective 
coupling field including mean-field effects is 



O n (r,i) = Q(r,i) + fs(r,t). 



(5) 



The collision integrals in Eq. (||) and Eq. (|J) are written 
explicitly in the Appendix. 

The center of mass and spin are coupled via the third 
and fourth terms in both Eq. (El) and Eq. (||). These 
terms involve spatial gradients over the density and scale 
like gn/ksT <C 1, and thus can be neglected. This allows 
us to make the further simplification that the center of 
mass and spin dynamics are decoupled. Since we are 
interested in the intrinsic collective modes of the system, 
we also assume that the field f2(r) = (0,0, A(r)) can 
be made to vanish after the spin wave is excited. This 
assumption is motivated by the JILA experiment p| , 
where the static inhomogeneous frequency splitting could 
be adjusted to zero after a short excitation time. The 
position dependence of A(r) determines the symmetry 
of collective mode excited. With these two assumptions, 
the kinetic equation for the spin distribution function 
then simplifies to 



da 

+ — • Vct - Vt/ oxt 
at m 



n - 9 Ft -. da 



(6) 



coll 



In this paper, we consider small amplitude spin oscil- 
lations around a fully-polarized state. It is convenient 
to define new spin coordinates (u',v',w') with w' being 
the direction of the initial spin polarization. In the ex- 
periment described in Ref. JlT| ], the spin oscillations oc- 
cur about the spin state polarized along the v axis. In 
this case one must make the mapping of the spin coor- 
dinates (u',v',w') = (w,u,v), corresponding to a cyclic 
permutation [ p3| . We then linearize the kinetic equation 
around the equilibrium state polarized along the w' di- 
rection (we drop the prime from our notation from here 
on): <7„o = ctdO = 0, and a w0 (r,p) = / (r, p), where the 
equilibrium distribution is 



/o(r,p) = exp<^ -j3 



2m 



+ U eKt (r) - no 



The equilibrium density given from Eq. (0) is 
n (r) = 



( 2 ^ 3 /o(r,p) = ^-exp[-/3C/ cxt (r)], 



(7) 



(8) 



where A t h = (2irh 2 /mk^T) 1 / 2 is the thermal de Broglie 
wavelength. We then substitute <?(r, p, t) = <?o(r, p) + 
<5<r(r, p, t) into (||) to obtain the linearized spin kinetic 



equation 



dSa 



+ — ■ VSa - VC/ cxt • V„(5<T 
m 

9 /Ft r ^ , v da 

- ~(S x da + SS x a ) = 



(9) 



The linearized form of the collision integral is discussed 
in the Appendix. 

The longitudinal spin dynamics is described by 



dSa v 



— ■ VSa n 



VC/cxt ■ Vp&T 

V 



dSa„ 



(10) 



coll 



Since the mean field term does not appear in (|10|) , collec- 
tive oscillations of the longitudinal spin only occur in the 
low-density collisionless regime in a trapped gas. In the 
high-density hydrodynamic regime, longitudinal modes 
become purely relaxational modes damped by the dif- 
fusion transport process. The crossover from the colli- 
sionless spin oscillation to the hydrodynamic relaxation 
mode in a longitudinal spin excitation was observed in a 
trapped 40 K fermi gas |24|, (although they were not 
excitations from the fully polarized state as considered 
here). In contrast, transverse spin waves behave collec- 
tively due to the mean field and will be the focus of this 
paper. For the transverse spin fluctuations, it is conve- 
nient to work with Sa± = 5a u ± iSa v . We then obtain 



dSa± 
dt 



m 



VSa± - VU C 



V p Sa± 
dSa± 



dt 



(ii) 



Although we are interested in a trapped Bose gas, it is 
useful to summarize earlier results on the theory of spin 
waves in a homogeneous gas |26| . In the long wavelength 
limit, where the gas can be treated in the hydrodynamic 
regime, the dispersion relation has the form M |26| 



-ik 2 v 2 f D 



(12) 



where v t h = y^ksT/m and td is a complex diffusive re- 
laxation time. For transverse spin oscillations, one finds 



i L , = (Tp 1 — ign/h), where the diffusive relaxation time 
is t d = [(32a 2 n/3)^7rfc B T/m] 



Due to the exchange 
mean field, the transverse spin behaves collectively, and 
in the limit where gnTn/Ti ^> 1, the dispersion relation 
has the form 



2;„2 



2m* 



to V gn J ' 



(13) 



where m* = (gn/2kBT)m is regarded as an effective 
mass. The k 2 dispersion relation is a universal result 
for ferromagnetic- like spin systems (T) . The longitudinal 
spin oscillations do not behave collectively, but rather ex- 
hibit a purely diffusive mode, with fry — > td in Eq. (J13) . 
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III. MOMENT METHOD FOR A TRAPPED 
SPIN- 1/2 GAS 

We now turn to spin waves in a Bose gas confined in 
the harmonic trap potential. Starting from Eq. (11), one 



can derive a general set of coupled moment equations 
associated with a set of polynomial functions Xi( r iP)'- 



dt 



m 



(V[/ ext ■ V p Xi 



+ M(s + Xi)o - (mxi)} = (xi)coii, 



where the moment variables are defined as 

{Xt) = Jj J dT J (2§)3 Xi(r,p)Sa+(r,p,t), 



(14) 



(15) 



where u>i is either u> z for an axial mode or u± for a radial 
mode. These are not a closed set of equations, since in 
general the mean field and collisional terms in Eq. j20| ) 
couple to higher moments. The hierarchy of moment 
equations can be truncated by assuming the explicit trun- 
cated form for the distribution Sa+: 



5ua 



h[oto + ol\x 1 + a 2 p Xi 



(21) 



The coefficients in the expansion can be related back 
to the set of moments using (|l5|): ao = 0, Oi\ = 
(mojf/k B T)(xi), and a 2 = {x2)/k B T. 

Using the explicit form for 5a + in Eq. ( pl| ) , the mean- 
field and collisional terms can be expressed in terms of 
the dipole moments. The resulting closed set of coupled 
moment equations, written in matrix form, is 



(Xi 



oil = I dr 



dp 



(2-n-ny 



dSa . 



dt 



(16) 



coll 



(Xi/0 = 



dr 



dp 



(2nhy 



- Xl (r,p)/o(r,p). (17) 



In general, moment equations are not closed since the 
mean field and collisional terms couple to higher mo- 
ments. Those higher moments can be truncated by ex- 
panding the fluctuations in the distribution function <5<7+ 
in powers of position and momentum. One can relate the 
coefficients in the expansion to the moments of the dis- 
tribution function (xi) to yield a closed set of equations 
that can be solved analytically. 

The choice of the functions Xi depends on the sym- 
metry of the spin-wave mode we are interested in. In 
the following subsections, we consider the dipole and 
quadrupole oscillations. 



A. Dipole mode 

We first consider the spin-wave collective mode with 
a dipole symmetry, which could be excited by a linear 
inhomogeneous frequency splitting, such as A(r) oc z. 
For our set of moments to describe this oscillation, we 
choose: 



Xi = Xi, X2 =p x Jm. 



(18) 



The subscript indicates the axis along which the excita- 
tion oscillates Xi € {x, y, z}. Within the moment method 
approximate treatment, the dipole modes along the three 
axis are completely decoupled. We also define xo = E 
related to the norm of 6a+; we set (xo) = 0, which is 
required from the conservation of the spin density. 

From Eq . (|14[) , the equations of motion for the mo- 
ments Eq. jlqjare given by 



S <Xi> = <*>. 
|(X2)+-?<Xi>-*f 



n X2) 



(19) 

(X2>coU, (20) 



dt x - Wdx ' 

where the coupling matrix Wd is 



W d = - 



-1 
2 (j D - jwmf) 



(22) 



(23) 



Here we have defined the vector of moments x = 
{(xi)j(X2)}- Three different frequencies appear in Wd'- 
the trap frequency u>i along the axis of oscillation, the 
mean field frequency lomf defined by the spatial average 
ujmf = J drgnl{v)/NTi 



WMF 



gn (o) 

2V2H ' 



(24) 



and the spatially averaged diffusion relaxation rate jo, 
the form of which is given in the discussion in the Ap- 
pendix. 

We now look for normal mode solutions x = Xo e ~ tut . 
Substituting this into Eq. ( p2| ) yields an eigenvalue equa- 
tion with two solutions. It is straightforward to show 
that the mode frequencies to obey the dispersion relation 



wf = 0, 



(25) 



where jo = Id — iwmf is the effective (complex) diffu- 
sion relaxation rate including the mean-field effect. The 
solution is given by 



1 

uj = — 
2 



-ijD ± \/ 4w, 



(26) 



We shall consider two limiting cases to obtain the scaling 
behavior of the frequency and damping of the modes. In 
the weak interaction, or collisionlcss, limit where tOi 3> 
|7d|, one has 



±u>. 



.Id 



(27) 



In the strong coupling, or hydrodynamic, limit where 
Wj <C |td|, one has 



-iujf/jD (low) 
-Z7D (high) 



(28) 
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These represent low and high frequency modes, with the 
low frequency mode having a higher Q ~ Rcw/Imw value. 
The low-frequency solution has the form of the diffusion 
relaxation rate with a complex diffusion coefficient. More 
explicitly, the dispersion relation in the strong-coupling 
limit takes the approximate form 



WMF 



MF 



. ID 
I 

WMF 



(29) 



The coefficients in the expansion can be related back to 



the set of moments using Eq. (|l5|): 



<4<X4> + (X3> + (x6>]/2fc B T, a x = (mw 2 /fc s T) 2 (xi)/2, 



a 2 



= mw 2 (x 2 }/(fc B T) 2 



03 



= ( X3 )/2(k B T) 2 



04 



(mwi/fc B T) 2 (x4)/2, a 5 = mw 2 (xsV^T) 2 , and a 6 = 
( X6 )/2(k B T) 2 . 

Following the same procedure as for the dipole mode, 
and obtain a closed set of coupled moment equations, 
written in matrix form as 



The scaling behavior of the real part of Eq. ( J29| ) can be 
recovered in a simple model based on the homogeneous 
gas result of Eq. . For low-frequency collective modes 
in a trapped gas along the Xi direction, the wave vector 
k is estimated as k ~ 1/-Ri, where Ri — mmi/k-QT is the 
size of the cloud along the Xi direction. 



B. Quadrupole mode 

In the JILA experiment fPH , a spin- wave collective 
mode with a quadrupole symmetry is excited due to an 
approximately quadratically varying frequency splitting 
A(r) oc z 2 . In principle, the oscillation may be excited 
along both the axial or radial directions, and so we take 
the following quantities for our set of moments: 

Xi = z 2 , X2 = zPz/m, X3 =P 2 z /m 2 , (30) 
X4 = r\, \5 = r_L • P±/m, xe = P 2 j_/m 2 , (31) 



where r± = x 2 + y 2 and p± = yP 2 + P 2 ■ We also 
define xo = 1, related to the norm of 5a, which we set to 
zero (xo) = 0. 

The six moment equations for the above quantities are 



dt X = W « X > 
where the coupling matrix W q in this case is 



(39) 



W a 



( 


-2 














\ 




1d/2 
2w 2 


-1 






































-2 
















-i 


ID 12 
2w 2 


-1 




V 





25^t 







/ 



(40) 



The vector of moments is x = {(xi)> (X2}, • • • ■> (X6)}- The 
tilde on the relaxation rates indicates the complex form 
7 = (7 — iwmf)- The quantities 77. and 7^, given explic- 
itly in the Appendix, are the spatially averaged axial and 
radial thermal relaxation rates, the difference of which 
^7 = 7t — 7t i s n °t zero m general. We note that W q 
is nearly block diagonal and that the axial and radial os- 
cillations are coupled only through collisions associated 
with Sjt- 

By substituting x = Xne~ liA)t into Eq. (|39|), we obtain 
an eigenvalue equation with six solutions; the dispersion 
law is determined from 



+ ^M-FaM = 0, 



(41) 



d(xi) 

dt 
d(X2) 

dt 

d(x 3 ) 
dt 

djXi) 

dt 
d(X5) 

dt 
d(xe) 

dt 



2(X2> - 0, (32) 
- (X3> + vlixi) - *V ("0X2) = (X2>coii, (33) 

+ 2uJ z(X2) + if [(S+X3)o - («0X3)] = (X3>coll, 



where 



2(X5 



0. 



(34) 
(35) 



- (Xe) + ^±(Xi) - ir( n oXs) = (X5>ooU, (36) 

+ 2wi(x 5 ) + «| [(5 , +Xe)o - (noXe)] = (Xe)coii, 

(37) 



Just as in the previous case of the dipole mode, we trun- 
cate the hierarchy by assuming an appropriate form for 
the distribution 



5a + = f [a + axz 2 + a 2 zp z + a 3 p 2 z 
+a 4 r 2 ± + a 5 r ± ■ p ± + aQp 2 ^]. 



(38) 



uj 2 - 4w 2 



2w 2 



\ 1 — iuj/-j T 

^(lo) = (w + i7r)( w + *7t) 



+ (42) 
(43) 



F 3 (u>) = 28^ ( w 2 - 2w 2 + iw-y 



( w 2 - 2uii + iw^Y ) . 



(44) 



In both the weak interaction (w, 3> |td|, IttI) anc l strong 
interaction (wj -C |td|j |7tI) limits, the axial and radial 
modes are uncoupled and are determined from T\{uj) — 
0, with T\{ui) given in (p^). In the weak interaction limit, 
one has three modes 



2uJ i 2 1^ lD +lT >' 



UJ 



- [2UO, 

WMF 



WMF \ l , 

- Z (7d 



77t- 



(45) 
(46) 
(47) 



FIG. 2: Dipole excitation. The surface shows how the spin 
component S v (z,t) varies with position and time. For visual 
clarity, we also show the contours projected below the surface. 



In the strong interaction limit, there is one (weakly 
damped) low-frequency mode and two (strongly damped) 
high-frequency modes. The low-frequency mode is given 
by u) = — Aiujf/jjj, which can be written approximately 
as 



UlM F 



2 / 
P V 



(48) 



This has the same scaling behavior with respect to the 
density, scattering length and temperature as the dipole 
mode result in Eq. (|29|). 



IV. NUMERICAL SOLUTION OF ID SPIN 
KINETIC EQUATION 

In this section we compare the predictions of the mo- 
ment method to a direct numerical solution of the spin 
kinetic equation. In Refs. |13, H, a one-dimensional 



model of the kinetic equation was presented and found to 
give very good agreement with experiments jl(J. A jus- 
tification that was given for the model is the separation 
of time scales between the axial and radial directions. 
In hindsight, the success of the one dimensional model 
can be further understood based on the moment method 
results in the previous section that the axial and radial 
modes are uncoupled in the linearized regime. 

We construct a one-dimensional model of the system 
by making the ansatz er(r,p,i) — a(z,p,t)ho(r±,p±) 
and then averaging over rj_ and p^. Here we take the 
static profile in the radial direction to be of Gaussian 
form h = ex.p[—(p 2 L /2m + mu!^r']_/2)/kBT]. We sub- 
stitute this ansatz into (Q) and integrate over the ra- 
dial phase space variables, which gives the following one- 



FIG. 3: Quadrupole excitation, as in Fig. 2 



dimensional model Boltzmann equation 

da p da dlf cxt da ^ _ da 
dt m dz dz dp dt 



(49) 



ID 



Here we have dropped terms that scale like gii/hsT. 
The collision integral in one dimension involves a phase 
space average in the radial direction da/dt\ 1D = 
■fxy ^/^Icoll / I X y wnere we have introduced the no- 
tation J xy ■■ ■ = J dr± J dp± ■ ■ ■ /(2irfc) 2 . The radial av- 
eraging introduces a scaling factor in the mean field 
terms, so that g — > g' = g/(2A 2 h ), where Ath is the ther- 
mal de Broglie wavelength, g' has the correct units of 
energy times distance required in our one dimensional 
model. 

Although the direct numerical simulation using the full 
expression for the one dimensional collision integral de- 
rived from Eq. ( p4j ) is technically feasible, we introduce 
a simple model for the relaxation 



da 

In 



1 



ID 



Tel (z) 



[a(z,p,t)-M(z,t)f Q (z,p)l (50) 



where t c \{z) = [16a1 2 n o( z ) V TtfosT / 'm] 1 is the radially 
averaged mean collision time, fo(z,p) = /o( r , p)/^0i an d 
M(z,t) — S(z,t)/nQ(z). Equation (|50j) contains the 
essential properties of collisions: (i) it vanishes when 
the distribution function has the local equilibrium form 
a(z,p,t) (x M(z,t)e~ p / 2mk ^ T ^ (ii) it conserves the spin 
density. We note that the form Eq. ( |50"| ) does not require 
the knowledge of the long-time equilibrium solution for 
S(r,t). 

We solve Eq. ( f49| ) numerically using a finite differ- 
ence, alternating-direction Crank-Nicholson routine. As 
a check on the numerics we monitor the integrated spin 
components J dzSi(z,t), which are conserved if = 0, 
and the integrated spin magnitude, which must satisfy 
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0.5 1 1.5 2 2.5 



density (10 13 cm" 3 ) 

FIG. 4: Frequency and damping rate of the dipole mode ver- 
sus peak total density no(0). The solid line is the low fre- 
quency solution obtained from Eq. ( p5| ) and the chain and 
dashed lines are the weak Eq. ( p7[ ) and strong Eq. (^9|) inter- 
action limits, respectively. The fill ed p oints are obtained from 
a direct numerical solution of Eq. mm for three different tem- 
peratures T: 600 nK (green circlej7800 nK (red diamond), 
and 1 mK (blue square) . 




0.5 1 1.5 2 2.5 

density (10 13 cm~ 3 ) 



FIG. 5: Frequency and damping rate of the quadrupole mode 
versus peak total density no(0). The solid line is the low 
frequency solution obtained from Eq. (^2j) (e g. T \ = 0) and 
the chain and dashed lines are the weak Eq. (W3) and strong 
Eq. (^) interaction limits, respectively. The filled points are 
obtained from a direct numerical solution of Eq. ( J49| ) for three 
different temperatures T: 600 nK (green circle), 800 nK (red 
diamond), and 1 mK (blue square). 



the relation 

| J dz{Sl + S 2 y + Si) = 2jdz J- (51) 

The total spin can decay to zero if there is an inhomogc- 
neous field tt(z,t) present (i.e. Vf2 ^ 0). 

To compare directly with the moment method results 
derived in the previous section, we take as an initial state 
a(z,t = 0) = a (z) + Sa(z), where a Q (z) = {0,0, fo(z)} 
and Sa(z) — {ReSa + (z), ImS<r + (z) 7 0}. We take 5a + 
of the form given in Eq. ( pl| ) for the dipole mode and 
Eq. ( |38| ) for the quadrupole mode. The coefficients oti 
are determined by diagonalizing the coupling matrices 
Eq. (^3|) and Eq. (^0|) respectively (in this section we 
consider only the low frequency excitations for a given 
symmetry) . 

In order to visualize the spatial form of the spin wave, 
we show the transverse spin component S v {z,t) in Fig- 
ures 2 and 3 as a function of position and time to show 
the symmetry of the dipole and quadrupole modes. We 
note the qualitative agreement between the contours of 
Figure 3 and the graph shown in Figure lb of Ref. [ pTj . 
We do not plot the S u (z,t) component, which has the 
same structure but is shifted 90° out of phase in time. 
The transverse spin {S u (z, t), S v (z, t)} at a given posi- 
tion z traces out a spiral that terminates at the origin as 



t — > co, whose overall diameter varies with position. The 
spin precession as a function of position for the actual 
experiment is shown in Figure 2 of Ref. |i~l| ] . 

We extract a frequency and damping rate from the 
numerical solution by calculating the dipole (z) or 
quadrupole (z 2 ) moment and fitting this quantity to a 
damped sine function of the form Aexp(— Bt) sin(Ci + 
D). The coefficients are obtained using a least squares 
fitting routine. The frequency to — C and damping rate 
7 = B are then compared to the predictions of the mo- 
ment method. 

In Figures 4 and 5 we plot the frequencies and damp- 
ing rates for the dipole and quadrupole spin waves. 
We take values for the physical quantities correspond- 
ing to the JILA experiment, where lu z /2tt = 7 Hz, 
and ojj_/2tt = 230 Hz. We find that the frequencies 
obtained from the numerical solution agree extremely 
well with the moment method prediction for the dipole 
mode, while the quadrupole mode shows only qualita- 
tive agreement. We also find that the frequencies in 
Fig. 4 and Fig. 5 are given as temperature-independent 
functions of the peak density, which is consistent with 
the moment results. In general, one can show that the 
streaming term of the linearized spin kinetic equation 
in Eq. (O) is invariant under the scale transformation 
T -> T',r -> (T'/T)^ 2 r,p -> (T'/T)^ 2 p for given 
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no(0), and thus the frequency of any given mode is tem- 
perature independent. In Figures 4 and 5 we do not show 
the damping as a function of temperature, which scales 
approximately as yT at fixed peak density, according to 
Eq. @. 

The most interesting feature occurs in the damping, 
where we see that the damping of the quadrupole mode 
is qualitatively different from that predicted by the mo- 
ment method. We postulate that this difference, which 
is largest at intermediate densities, is due to Landau- 
type damping arising from the mean-field coupling of the 
collective mode to higher excitations. Our moment ap- 
proach using the simple truncated form for the distribu- 
tion function cr + , as given by Eq. (^Tj) or Eq. (|38|), does 
not account for this effect |23|| . 

For a homogeneous Bose gas, Oktel and Levitov 
worked out the spectrum of spin waves in the collisionless 
regime pl| using a linear-response theory with the ran- 
dom phase approximation. They find that the mean-field 
coupling gives rise to Landau-type damping, which is in 
addition to the collisional damping due to spin diffusion. 
The damping rate is given by 



1L 



V^g 2 n 2 _ ( gn /hkv th ) 



h 2 kv 



(52) 



th 



A rough estimate of Landau damping in a trapped gas 
may be given by using k ~ l/i?th in this uniform gas 
result. However, this simple estimate predicts a damp- 
ing rate that is about an order of magnitude larger than 
that shown in Figure 5. Apparently one needs to work 
out the linear response calculation taking into account 
the excitation spectrum explicitly in a trap potential in 
order to obtain a better quantitative model for Landau 
damping. 

In Ref. we compare the numerical solution pre- 
sented in this paper to experimental data, and find ex- 
cellent agreement for both frequency and damping. It 
is important to realize that in our numerical calculation, 
we have used a relaxation time approximation to treat 
the collision integral. In this model, one has some free- 
dom to choose an appropriate relaxation time; there are 
a few different candidates, but here we have used the spa- 
tially dependent form of the mean collision time given in 
Eq. (50). Although this is a reasonable choice, it was not 
a priori obvious that it would result in the best agree- 
ment with experimental data. 

We also remark that our results are strictly valid for the 
case where the inhomogeneity A(r) is vanishingly small 
(though we have not carried out the full linear response 
theory, we envision A(r) as playing the role of an exter- 
nal perturbation used to excite the spin wave). In the 
JILA experiment [jll), the effect of applying a constant 
static A(r) during the entire spin wave oscillation was 
investigated as a function of the magnitude of the per- 
turbation. For a large inhomogeneity, the response of the 
system is nonlinear and the mode frequency is modified; 
the earlier JILA experiment JlO| seems to reside in this 
regime. 



V. CONCLUSION 

In this paper we have studied spin waves in a dilute 
noncondensed Bose gas of two level atoms. Our main 
new contribution is that we have treated an inhomoge- 
neous system held in a harmonic trap in order to de- 
scribe related experiments on spin waves JlT| | . We ap- 
plied the moment method technique for dilute trapped 
gases, which lead to closed form solutions for the frequen- 
cies and damping of dipole and quadrupole modes. As a 
test of the validity of our moment model, we compared 
the results to the numerical calculation of a one dimen- 
sional model of the spin Boltzmann equation and found 
very good agreement overall. The main discrepancy be- 
tween the moment model and the numerical calculation 
occurred in the damping of the quadrupole mode. We at- 
tribute this to the Landau damping, which is contained 
in the numerical approach, but absent in the moment 
model. 

Although the salient features of spin waves in a non- 
condensed dilute gas seem to be well described by our 
theory, which shows excellent agreement with experimen- 
tal data presented in Ref. , the situation for this sys- 
tem just below the Bose-Einstein condensation temper- 
ature [ fl8| remains largely unexplored. Many zero tem- 
perature properties of a spin-1/2 condensate have been 
studied, when the thermal component is absent, such as 
spin winding (^?J and vortex-like spin textures [^8[ |29| . It 
should be noted, however, that at zero temperature the 
condensate does not support spin waves of the type con- 
sidered in this paper, since the exchange mechanism does 
not occur in the condensate. It is very interesting that 
the spin-1/2 thermal gas exhibits strong collective behav- 
ior in the spin dynamics due to the exchange effect, even 
in the collisionless regime; this is in sharp contrast to the 
behavior of density fluctuations of a single component 
thermal gas, where the mean field of the noncondensate 
plays a very minor role (mainl y a s a source of damping 
of the condensate excitations) |30[ [II) . At finite temper- 
atures, the spins of the condensate and thermal gas will 
interact strongly. For future studies, it will be interesting 
to investigate how long-lived spin textures in the conden- 
sate are affected by the thermal gas, which in principle 
can itself support (probably short-lived) spin textures. It 
will also be important to understand how the condensate 
modifies spin waves in the thermal cloud ftn\ . 
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APPENDIX: COLLISION INTEGRALS AND 
RELAXATION TIMES 

The collision integrals appearing in the kinetic equa- 
tions Eq. (^) and Eq. (Q) for the center of mass and spin 
distribution functions are given by 



df 

at 



coll 



7T5 

h 



dp 2 



dp 3 



(2irh) 3 J (2tt?i) 3 
5(e p + e P2 + e P3 + e P4 

)5(P 

X{3[/( P3 )/(P4)-/(P)/(P2)] 

+a(p 3 ) • ct(p 4 ) - ct(p) • a(p 2 )}. 



P2 



P3 



da 
~dt 



h 

5(e p 



dp2 



dp 3 



dp 4 



(2irh) 3 J (2irh) 3 

vi ~ e P3 - e P4)^(P + P2 
x{3/(p 3 ) ( ?(p4)+a(p3)/(p4) 
-/(pMP2)-3a(p)/(p 2 )}, 



P3 



P4) 

(53) 

P4) 

(54) 



where e p = p 2 /2m. Here we neglect a principal value 
contribution, which gives a second-order correction to 
the free streaming evolution, and we take all scattering 
lengths <2y to be equal - a reasonable approximation for 
87 Rb. This approximation results in the conservation of 
spin density during collisions, i.e. J dpdS / dt\ co \\ = 0. 
When the small differences in scattering lengths are ac- 
counted for, the transverse spin decays slowly. For 87 Rb, 
this contribution to the "T2" lifetime is of the order of 

10 Sg|. 

For the linearized form of the spin collision term ap- 
pearing in Eq. (^), it is convenient to express the fluctu- 
ation of the spin distribution function as 



Sa(r, p, t) = f {r, p)ip(r, p, t). 



(55) 



Then the linearized collision integral is written as 
da/dt\ coll ee L$\ 



L bP] 



Kg 
h 

x 6(e p + 



dp2 



(2-Kh) s 



dp 3 



dp4 



l>2 



C P3 



(2tt^) 3 
ep4)^(P + P2-P3-P4) 



X /0(P3)/ (P4){2[V(P 4 ) - tf(p)] 

+ [^(P 4 )+^(P 3 )-^(P2)+^(P)]}- 



(56) 



Using the approximate forms of the distribution func- 
tion Eq. (|Tj) and Eq. ( |38"| ) in the linear collision operator 
and taking moments, we find that the collisional contri- 
butions to the moment equations are given by 



(Pxi)coll = -iDiPxi), 



{zpz) co\\ = 



(rj. ■ P±) 



coll 



Id 



r± ■ p±), 



(57) 
(58) 

(59) 



(^)coii = -1t(p 2 z) - $1t(p 2 ±), 



(pi)coll = -7T(pi)-2^7T^). 



(60) 
(61) 



The various relaxation rates are given by spatial average 
of the following spin transport relaxation times: 



"o(r) . 
r D (r) ' 

n {r) = 
r T (r) " 

n Q (r) _ 
T^(r) 



1 



dp 



3mfc B T J {2Trhf 

dp 1 

{2irhf 2{mk B T) 2 

dp 1 



pL[p], 



(62) 



P 2 z L[pI}, (63) 



{2irh) 3 4(mfc B T) : 



pirn 



(64) 



The associated spatially averaged relaxation rates are 
given by 



Id 



It 



1 

N 



—— / dr 



n (r) 



1 



1 



1 

N 



dr 



It 



1 

N 



= — dr 



r D {r) 


2V2~t d {0) 


n {r) 


1 


r T (r) 


2V2r|(0)' 


n (r) 


1 



T T (r) 2\/2r^(0)' 



(65) 



(66) 



(67) 



For detailed calculations of various transport relax- 
ation times in a trapped Bose-condensed gas, we refer 
to Ref. |S2| . It is straightforward to generalize those cal- 
culations to work out the spin transport relaxation times 
defined above. We find 



4 \r) 



3 

5 ' 



(68) 
(69) 
(70) 



where r c i(r) = [32no(r)a 2 (7rfcBr/m) 1 / 2 ]~ 1 is the mean- 
collision time. The spatially averaged relaxation rates 
are given by 

1 z 3 JL 7 2 

Id = g7oi 5 7t = g7ci, 7t = 7g7ci, 07t = _ Tg7ci, 

(71) 

where 7 c i is the spatially averaged mean collision rate 

1/2 



16 fn \ 2 ( nk BT 
7cl = -j=n (Q)a 



(72) 



In order to compare the moment method directly with 
our numerical solution of the ID kinetic equation, we also 
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worked out the three relaxation rates within the simple 
relaxation time approximation 



djf 
dt 



a(v,p,t) - M(r, t)/ (r,p) 



(73) 



coll T d( r ) 

This simple formula leads to the same collisional con- 



tributions in the moment equations as given from the 
original collision integral, with all the transport relax- 
ation times being replaced with r c i(r), so that one has 
Id = It — It — 7ci- Thus, within the relaxation time 
approximation, oscillations in the axial and radial direc- 
tions are completely uncoupled since Sjt = 0. 
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